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The updated CLEO measurement of the electromagnetic penguin process B — » 
X s 7 is presented and compared to the ALEPH result. Implications on new physics 
are discussed and the status of recent searches for B — > X s l + l~ modes is given. 



1 Introduction 

Flavor-changing neutral currents (FCNC) are forbidden in the Standard Model 
(SM) at tree level. They are, however, induced at higher orders via penguin 
processes or box diagrams. One such process is B — > X s "f, which is mediated 
by electromagnetic (em) penguin loops. The lowest-order SM diagrams are 
shown in Figure [l| The dominant contribution arises from the magnetic dipole 
operator O7. QCD corrections, however, cause mixing with the other oper- 
ators 0\ to 06 and 0%. The dominant short-distance effects are separated 
from the non-perturbative long-distance parts into perturbatively-calculable 
renormalization-scale-dependent Wilson coefficients. In B — > X s j all pertur- 
bative contributions can be absorbed into a single effective function, (/u), 



which reduces to cf^^ (5.0) = -0.300 in leading-order QCDlil However, a 
25% uncertainty in the decay rate arises from varying the renormalization scale 
between and 2m;,. In next-to-leading order (NLO) calculated recentlyB 

this uncertainty is reduced by a factor 1 .5-213 nn 
In the Standard Model the inclusive B — > X s j decay rate is aa 

T(B - X sl ) = ^G%a em I V tb V* | 2 (| C e 7 ff (^) | 2 +A), (I) 

where nib, Gf, ot em , A, and VtbV t * denote the &-quark mass, the Fermi con- 
stant, the em coupling constant, the bremsstrahlung contribution and Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements for t — > b and t —> s couplings, 
respectively. The B — > X s j branching fraction is calculated by scaling the 
semileptonic decay rate, yielding B(B -> X s ~/) = (3.28 ± 0.33) x I0" 4 in NLO 
QCDcl Since the branching fraction is sensitive to additional contributions from 
e.g. a charged Higgs boson, a new gauge boson or a charged supersymmetric 
(SUSY) particle in the loop, B — > X s ^ is well-suited to probe new physics. 
A branching fraction measurement e.g. can be used to constrain the SUSY 
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Figure 1: Lowest-order Feynman diagrams for B — ► X s ~f in the SM. 



parameter space. CLEO was tha first experiment to observe this process both 
in the exclusive K*j final stated and the inclusive channelQ 

2 The recent CLEO B -> X s -f Measurement 

The inclusive B — > X s j analysis has been updated recently by CLEO using a 
sample of 3.3 x 10 6 BB pairs. The photon energy, which is nearly half the B 
mass (< E 1 >~ 2.4 GeV), is the highest observed in any B decay. The Fermi 
motion of the &-quark inside the B mesons and the B meson momentum in 
the lab frame yield a Doppler-broadened photon line. According to spectator 
model predictions, 85.% - 94% of the photons are in the 2.1 GeV < E 1 < 
2.7 GeV energy range! While backgrounds from B decays are small, those from 
continuum processes, such as e + e~ — > qq~/ and e + e _ — > qq, are significant. To 
minimize these backgrounds an event shape analysis is combined with exclusive 
final-state reconstruction, considering only high-energy photons as candidates 
that are inconsistent with originating from a 7r° or r/ decay. 

At the T(45) the event shape of hadronic BB decays is spherical, while 
that of qq continuum is back-to-back. In B — > X s j events the spherical sym- 
metry is distorted by the hard photon. The event shape resembles more that of 
qq continuum with initial-state radiation (ISR). The ratio of second-to-zeroth 
Fox- Wolfram moments R2 and the sphericity S± are used to separate the sig- 
nal from qq continuum. To suppress ISR continuum background 2 variables 
are evaluated in the rest frame of the e + e~ following the radiation of a hard 
photon: the ratio of Fox Wolfram moments R' 2 , where the photon is excluded, 
and cos 9', where 9' is the angle between the photon and the thrust axis of 
the rest of the event. In addition, the energy flow inside 20° and 30° cones 
parallel and antiparallel with the photon direction are calculated. With help 
of a neural network the 8 variables are combined into a single shape variable 
r s h, which favors values of r 8 h ~ +1 for B — > X s j signal and r s h ~ —1 for 
continuum background. 
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For additional suppression of continuum background the X s system is re- 
constructed in exclusive final states. Considered are modes containing either 
a K^ or a K s — > n + ir~ plus 1-4 7r's with at most one n° in addition to the 
hard photon. A global x 2 based on the beam-constrained mass Mb, the energy 
difference AE, the dE/dx of the charged particles as well as the Kg and ir° 
masses is used for candidate selection. For the combination with the lowest \ 2 
a second shape variable r c is constructed using a neural network. The input 
variables are r s h , the x% va l ue obtained from Mb and AE alone and | cos Ott | , 
where Ott is the angle between the thrust axes of the candidate B and the rest 
of the event. Again signal events peak at r c sa +1, while continuum back- 
ground peaks at r c ss — 1. Figures §a,c show the r c and r s h distributions for 
signal and qq continuum obtained from Monte Carlo (MC). 




Figure 2: a) Distributions in the shape variable r c for MC samples of B — > X s -y and qq 
continuum for events with x% < 20; b) event weighting function for r c ; c) corresponding dis- 
tributions in shape variable r s ^ for events failing reconstruction; d) event weighting function 
for r sh . 

To optimize the signal-to-background ratio, the candidates are weighted 
by to = s/ [s + (1 + a)b] , where s and b are the expected signal and background 
yields for a given value of r c and a is the luminosity scale factor between 
the T(45) and the off-resonance data samples. Since just ~ 20% of signal 
events are exclusively reconstructed, a weighting based on r s h alone is used for 
selected events which fail the exclusive reconstruction (\ B > 20). The event 
weighting functions are displayed in Figures ||b,d. 

Backgrounds from BB decays, where the photon originates from a 7r° or an 



3 



r\ decay, are estimated from Monte Carlo. However, the tt° and 77 momentum 
spectra are adjusted to those observed in the data. The spectra are obtained 
by treating the ir° or ij like a photon, performing the same analysis as discussed 
above for photons using weights w. This procedure helps to account for inac- 
curacies in the b — > cW~ and b uW~ event generation and for omissions of 
modes such as B — > X s g. Further analysis details are given in reference [9]. 

Figure ||a shows the weighted photon spectrum and the estimated contri- 
butions from BB and continuum backgrounds. The background-subtracted 
photon spectrum plotted in Figure ||b shows a clear signal in the 2.1-2.7 GeV 
region. The observed shape is consistent with the spectator model prediction!] 
The weighted event yields in two bins below and one bin above the signal 
region are consistent with the background parametrization. The background- 
subtracted weighted event yield in the 2.1-2.7 GeV window amounts to 92.21 ± 
10.26 s tat ± 6A6 sy s events. The systematic error consists of a ±1% uncertainty 
in the qq continuum subtraction and a 20% uncertainty in the BB subtraction, 
both added in quadrature. 

The weighted efficiency in the 2.1-2.7 GeV is determined from the Ali- 
Greub modeip which includes gluon bremsstrahlung and higher-order radiative 
effects. The Fermi momentum of the &-quark and the mass of the spectator 
quark are varied such that on average < mt, >= 4.88 ± 0.1 GeV. For the 
hadronization process both the JETSET algorithm H and a mix of K* reso- 
nances_are used such that the X s mass distribution predicted in the Ali-Greub 
model is reproduced. A weighted efficiency of e w = (4.70 ± 0.23 ± 0.23 ± 
0.04 ± 0.33) x 10 -2 is determined in the signal region. The 4 errors represent 
uncertainties due to spectator model inputs, hadronization of the X s system, 
the production ratio of B + B~ and B°B° and detector modelling, respectively. 

Combining the background-subtracted weighted event yield with the 
weighted efficiency results in a branching fraction of 

B(B -» X sl ) = (3.15 ± 0.35 stQt ± 0.32 sys ± 0.26„ lO£ie; ) x 10~ 4 . (2) 

This measurement is in good agreement with the Standard Model prediction. It 
supercedes the previous CLEO result, since 60% additional data are included 
and improved analysis techniques are used. The 95% confidence level (CL) 
limits including all systematic effects are 

2.0 x 10~ 4 < B(B -> X sl ) < 4.5 x 10~ 4 . (3) 
3 The B -> X sl Result from ALEPH 

At LEP the b and b quark are separated into opposite hemispheres. The 
B — > V s 7 hemisphere has low multiplicity, a single displaced vertex and a high- 



4 




Figure 3: a) Observed weighted photon energy spectrum and background contributions from 
CLEO; b) background-subtracted photon yield and predicted shape of B — » X s ~i signal. 

Figure 4: top) Observed photon energy spectrum and background contributions from 
ALEPH; bottom) background-subtracted photon yield from multivariate fits. 



energy photon, whereas the opposite hemisphere contains a typical b hadron 
decay that is used for a 6-tagging. In a sample of 8.9 x 10 5 66 pairs, ALEPH 
selects B — * X s j candidates with an inclusive reconstruction algorithm. First, 
probabilities P to originate from a B — > X s j decay arc defined for the candi- 
date photon and the other particles (charged tracks, 7r°'s, Kg's and K®'s) in 
the candidate jet, using rapidity, momentum and impact parameter. Second, 
the candidate photon and up to 8 particles in decreasing order of probability 
P are added until an invariant mass of the system closest to the B mass is 
obtained. In addition to standard b hadron selection criteria, the candidate 
photon has to be inconsistent with a tt° decay; the 77 invariant mass when 
combined with another photon has to lie above 0.2 GeV. The photon energy 
has to exceed 10 GeV and 9*, the angle between the photon in the rest frame 
of the jet and the jet direction has to satisfy cos 9* < 0.55. While the signal is 
uniformly distributed, the background peaks at cos 9* = +1. The mass of the 
candidate B has to be within 0.7 GeV of the nominal B mass and the mass of 
the X s system has to be below 4 GeV. 

The resulting sample is devided into 8 subsamples, using the major axis of 
the electromagnetic shower ellipse 07 to separate merged 7r°'s from photons, the 
jet energy Ej et , and the probability P^em that all charged tracks in the opposite 
hemisphere are consistent with originating from the primary vertex. This 
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procedure provides an effective separation of signal photons from backgrounds 
that originate from 7r°'s, final state radiation (FSR) and other sources. The 
signal is enhanced in the subsample with Ej et > 32 GeV, 07 < 2.3 cm and 
— logP^^ > 2.2. The photon energy E* in the jet rest frame is plotted 
for this subsample in Figure [|. An excess is seen in the 2.2-2.8 GeV region. 
A binned log-likelihood fit is performed in all eight subsamples to extract the 
number of signal events and background events due to FSR, 7r°'s from b decays, 
7r°'s from non-6 decays and other sources. Of the 1560 candidates 69.4 ± 19.7 
are found to be signal events. The largest backgrounds result from ir° decays. 
The efficiency for selecting B -^-X s ~f signal is e = 12.8 ± 0.3%. Further details 
are discussed in the publicationt2l 

From the yield in the 2.2-2.8 GeV energy interval one obtains an inclusive 
branching fraction of 

B(B -> X al ) = (3.11 ± 0.8 ± 0.72) x 10~ 4 . (4) 

This is in good agreement with the CLEO result and the SM prediction. The 
largest systematic error contributions are due to background shapes (0.462), 
background MC statistics (0.376), background MC composition (0.185), and 
energy calibration uncertainties (0.182). 

4 Implications on New Physics 

New physics processes can enhance the B — » X s j decay rate. For example, in 
SUSY models one expects contributions from a charged Higgs and charginos in 
the loop as shown in Figure |. These yield new contributions C™ ew (Mw) and 
Cg ew (M\y) to the SM Wilson coefficients at the Af^-scale. Due to operator 
mixing, B — » X s j then limits the possible values for C" ew (Mw)- To study 
the effect of B — > X s ^ on SUSY, solutions in the SUSY parameter space have 
been generated using the minimal supergravity model (SUGRA)Ej The ranges 
of the input parameters are < to < 500 GeV, 50 < nfyi < 250 GeV, 
-3 < A /m < 3, 2 < tan/3 < 50 and m t = 175 GeVB Each solution 
is only kept if it is not in violation with present SLC/LEP constraints and 
Tevatron direct sparticle production limits. For these solutions the ratios Rj — 
C? ew (M w )/C$ M (M w ) and R s = C£ ett '(M w ) / '§ M ' (M w ) are calculated. The 
resulting scatter plot in the R7 — Rs plane is shown in Figure ||a. Many 
solutions are already in conflict with the present B — > X s j measurement from 
CLEO (region inside solid bands). 

To demonstrate the reach of B — > X s j in probing SUSY in comparison to 
that of high-energy colliders, 5 points in the. SUGRA model chosen at Snow- 
mass 1996 for an NLC study are examinedl!3 Point 3 is the common point for 
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Figure 5: Feynman diagrams for new physics contribution to B —* X s -f. 




Figure 6: Scatter plot of Kg vs i?7 for a) solutions in the SUGRA model and b) for 5 NLC 
SUGRA points. The allowed region from the CLEO measurement lies inside the 2 sets of 
solid diagonal bands. The dashed lines indicate a potential 10% measurement. 



a comparison of SUSY studies at NLC, LHC and upgraded Tevatron. The 
sparticle mass spectra for these 5 points are obtained by SUGRA relations and 
their contribution to-,£> — > X s ^f can be computed. Figure ^> shows the results 
in the Rj-R% planets Also shown are the constraints obtained from fits to the 
present CLEO data (solid line) and a future 10% measurement assuming a 
SM value for the branching fraction (dashed line). While point 4 is already 
excluded by the present CLEO measurement, all but point 5 lie outside the 
dashed diagonal band. Thus, rare B decays provide a complementary approach 
to high-energy colliders in searching for SUSY. 
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Figure 7: Lowest order Feynman diagrams for B — ► X s l + l in SM. 



5 Study of B -> X^+Z decays at CLEO 

Another interesting and related process is B — > A S Z + Z _ , where the photon 
is replaced with a lepton pair. The lowest-order SM diagrams are shown in 
Figure ^. In addition to the magnetic penguin operator O7, the electroweak op- 
erators Og and Oio contribute, yielding 3 scale-dependent Wilson coefficients: 
Cj(^i), Cg(fi) and Cio(n). Due to QCD effects, which are smaller here than 
in B — » X s j, again operator mixing occurs. Their contributions are absorbed 
into the effective Wilson coefficients Cy e ^(/z) and Cg^(fi). The decay rate 
has been calculated in NLO. However, to maintain a scheme-independent re- 
sult, only Cg(fi) is calculated in NLO while LO is retained for all other Wilson 
coefficients. New physics can affect any of these 3 Wilson coefficients. Due 
to the additional a em coupling, the branching fractions predicted in SM are 
almost 2 orders of magnitude smaller than those for corresponding B — > X s ~f 
modes. 

Using a sample of 3.3 x 10 6 BB pairs CLEO has searched for both in- 
clusive X s l + l~ channels and 8 exclusive final states, where a if or K* recoils 
against or e + e~. A detailed discussion of the inclusive analysis is given 

in reference [16]. In the exclusive analysis a Fisher discriminant method is 
used in addition to standard B selection criteria to suppress continuum back- 
ground. It is constructed from the distributions of thrust and sphericity of 
the B candidate, the B direction in the lab frame and the angle between the 
thrust axes of the candidate B and the rest of the event. The candidate B 
is selected by a x 2 method, which includes AE, the masses of intermediate 
resonances, and both the dE/dx and time-of-flight information of all charged 
tracks in the final state. The combination with the lowest \ 2 is retained. Since 
double semileptonic decays are a dominant background source, the missing en- 
ergy of the event is constrained to < 20%. To remove contributions from the 
color-suppressed B — > modes, events with dilepton masses near the tp 

and ip' mass peaks are excluded. Signal regions are defined in the Mb - AE 
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Figure 8: Branching fraction upper limits of B — > X s l + l~ modes at 90% CL from CLEO 
(solid points) in comparison to a range of SM predictions (solid lines), recent results from 
CDF (squares) and previous upper limits (triangles). The numbers on the right-hand side 
show the present CLEO limits in units of 10~ 5 . 



plane by ±2. 5(7 a/ x ±2.5ctae wide boxes, where <jm, cab denote Mb and AE 
resolutions, respectively. Regions outside 3.5ctm and a^E are used for back- 
ground determination. No signal is observed in any of the 8 final states. The 
observed yields are consistent with expected backgrounds. The total efficiency 
is 4 — 5% for modes containing a Kg or a K* + , 8 — 10% for modes with a K*° 
and 17 — 23% for modes with a K + . Further details are found in reference [17]. 

The resulting 90% CL upper limits are platted in Figure || together with-. 
those of previous studies,ll3 recent CDF resultsEj and a range of SM predictionsp 
The K*°l + l~ upper limit has been set for TO;+;- > 0.5 GeV, to reduce the 
contribution from . The new CDF and CLEO upper limits provide sig- 

nificant improvements, although they still lie above the SM predictions. The 
CDF upper limit closest to the SM prediction is that for K lying just a 
factor of 1.2 above. The most sensitive upper limit from CLEO is that for the 
combined K*°l + l~ modes, which lies a factor of 2.2 above the SM predictions. 



9 



6 Conclusion 



The present B — > X s 7 measurement is in good agreement with the NLO SM 
prediction. The combined statistical and systematic errors are 15%, while the 
model dependence is 8%. Even with the present accuracy, B — ► X s -f already 
provides important contraints on the SUSY parameter space. With the start 
of the BABAR, Belle and CLEO III next year, the accuracy of the B — > X s j 
branching fraction measurement will be significantly improved shortly. Suffi- 
cient B mesons will be produced in these experiments to make an observation of 
B — ► X s l + l~ modes possible. Thus, em penguin decays may actually provide 
the first hints on physics beyond the SM. 
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